Measurement-driven quantum evolution from a known state 
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We study the problem of driving a known initial quantum state onto a known pure state without 
using a unitary evolution. This task can be achieved by means of von Neumann measurement 
processes, introducing A'' observables which are consecutively measured in order to approach the 
state of the system to the target state. We proved that the probability of projecting onto the 
target state can be increased meaningfully by adding suitable observables to the process, that is, it 
converges to 1 when N increases. We also discuss a physical implementation of this scheme. 

PACS numbers: 03.67.-a, 03.65.-w 



I. INTRODUCTION 

The problem of controlling quantum systems has been 
a renewed subject of study. Quantum computing is 
basesed on the existence of a set of universal quantum 
gates which, concatenated, allow one to implement any 
unitary transformation within a fixed level of accuracy. 
These quantum gates are implemented through the con- 
trolled manipulation of the interactions among different 
physical systems. Quantum communication protocols, 
such as quantum teleportation entanglement swap- 
ping PI and dense coding, also require the precise ap- 
plication of some unitary transformations in a finite set 
of transformations. A related problem has also been ad- 
dressed in the contex of quantum control ^1§\. There, the 
goal is to drive the evolution of an initial, possibly mixed, 
state to a state having a predetermined expectation value 
of some observable. This evolution is also considered to 
be unitary. 

In this article we study the control of quantum sys- 
tems in the case where it is not possible to resort to uni- 
tary transformations. Our main goal is to map a known 
quantum state onto another known state via a sequence 
of measurements with the highest possible success prob- 
ability, that is, a controlled evolution via measurements 
only. It has been shown |l9l | that the mapping of an 
unknown quantum state onto a known pure state can 
be optimally implemented with the help of two observ- 
ables only. In this case, the highest success probability 
is achieved when the eigenstates of the two observables 
define mutually unbiased bases. It has also been shown 
that, when the system subjected to the measurements is 
affected by a decoherence mechanism, only one observ- 
able is required [20|. 

Here we study the case of driving by von Neumann pro- 
cesses [21] a known initial state makeing use of more than 
two observables. First we analyze the problem of two ob- 
servables. Thereafter, we show that a new observable can 
be added in order to achieve a further increase in the suc- 
cess probability. By means of numerical simulations we 
show that the success probability rapidly approaches the 



unity when the number of observables increases. 



II. DRIVING THE EVOLUTION BY TWO 
OBSERVABLES 

Let us start by supposing that a quantum system is in 
a known p state. Our goal is to drive the system to the 
known \C,) target state by measurements only. If we mea- 
sure the C observable, whose eigenstates are {|C): IC±)}i 
the probability of projecting to the |C) target state is 
Pd — (CIpIO- Natural questions arise: is it possible to 
increase this direct probability pd by making use of an in- 
termediate measurement of another observable 91 And, 
if it is possible, then how is the relation among p, C, and 
which maximizes such probability? 

So, in order to approach the state of the system to 
the IC) target state |22, we first measure an observable 
which has the {|0i), |li)} eigenstates. As a second step 
we perform a measurement of C,. Thus, the probability 
of reaching the |C) target through one eigenstate of 6 
followed by a measurement of C is given by 

PM = (0iH0i)|(0i|OP + (iiHii)l(iilC)p. (1) 

Making use of the normalization of p and and of the 
orthonomalization of |0i) and |li) the previous expres- 
sion can be cast in the form 

p,.s = (CIpIC) - 2 KiilcXClOi)!' ((CIpIC) - (a|pia)) 

+ (2|(0i|C>P - 1) [(Oi|C>(a|Oi)(C|p|C±> +C.C.] . (2) 

The second term at the r.h.s. of Eq. (|2Jl contributes to 
increase pi^s with respect to pd when (Cj-|p|C±) is higher 
than (CIpIC)j otherwise it helps to decrease pi^s with re- 
spect to Pd. Meanwhile, the third term plays a role when 
the p initial state has non-diagonal elements different 
from zero in the C representation. If p = //2, being / 
the identity, then pis = 1/2, so that pi,s is independent 
of the choice of 6*. If the p initial state is diagonal in the 
C representation, then, when (C|p|C) < (Cj^IpICj-); it re- 
quires a 6 observable unbiased to Q in order to optimize 
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FIG. 1: Probability of success pi s as a function of |(Oi|C)|'^ 
with (CIpIO = (solid), (CIpIC) = 0.5 (dash), and (ClpjC) = 
0.9 (dot). Horizontal lines are the respective pd- 



the process, whereas when (C|p|C) > (C±I/o|Cj_)i there is 
not any 9 observable which allows to increase pi^s over 
the value of pd- 

The third term of the r.h.s. of Eq. Q contributes 
maximally to pi,^ when arg((C|/9|C±)(C-L|0i)(0i|C)) = 
and KOilOP > 1/2, or arg((C|p|a)(a|Oi)(Oi|C)) = ^ 
and |(Oi|C)p < 1/2. Since both cases are symmetric with 
respect to |(Oi|C)P = 1/2, in the following we consider 
only the latter. 

FigureQlshows pi^s as a function of |(Oi|C)P for differ- 
ent initial values, say: (CIpIC) = (solid Hue), (CIpIC) = 
0.5 (dashed line), and (C|p|C) = 0.9 (dotted hue). In all 
of these cases we have consi dered complete in itial coher- 
ence, this is, |(C|p|a)l = V{C\p\0{C±\p\C±)- The hori- 
zontal lines are the respective pd = (CIpIO- We can see 
that for the considered initial conditions there is an in- 
terval of |(Oi|C)P where pi^s is higher than its associated 
Pd, and there is a particular value of |(Oi|C)p for which 
Pi^s is maximum. 

Let us examine what happen for a more general ini- 
tial condition, i.e., a p initial state with < |(C|p|C±)l 1^ 
v/(C|p|C>(a|p|a>, where (C|p|C) ^ 1- We look for condi- 
tions under which the probability of success of the mea- 
surement process M((^)M{9) is maximum. Optimizing 
Eq. with respect to |(Oi|C)P, one finds that the max- 
imum value Pmax of pi^s is 



(CIpIC) 



4(1 + ^): 



(3) 



with 



R = VO^lWX\p\C)^W+l^, (4) 
where we have defined the 7 coefficient by the equality 



mp\a)\ = W{c\p\o{a\p\a), o<7<i, 

Fig. (j21a) shows the maximum probability JSJ as a func- 
tion of (CIpIC) for different values of 7: 7=1 (dot-dashed 
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FIG. 2: (a) maximum probability of success pmax as a func- 
tion of (CIpIC) for different 7 values: 7 = 1 (dot-dash), 7 = 0.7 
(dot), 7 = 0.4 (dash), and 7 = (solid), (b) |(Oi|C)P com- 
ponent as a function of (C|p|C) for different 7 values: 7 = 1 
(dot-dash), 7 = 0.7 (dot), 7 = 0.4 (dash), and 7 = (solid). 



line), 7 = 0.7 (dotted hue), 7 = 0.4 (dashed hue), and 
7 = (solid line) . The diagonal solid line corresponds to 
Pd- Notice that for all 7 ^ the optimal probability 

Pmax 

exceeds pd for all values of (CIpIC)- Also, larger values of 
7 result in larger values of 

Pmax • 

The |0i) eigenstate of the 9 observable which optimizes 
Pi^s has a component on the |C) target state given by 
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1(0.101^^^ 1-^Vi 



2(C|p|C)-i 

R 



(5) 



Fig. ((21b) shows the square module Q of the |0i) state 
component onto the target state |C) as a function of the 
initial probability of the |C) state for different 7 values. 

When the initial state is pure (7 = 1), p = |'0)(V'|j the 
square module of the |0i) state component on the target 
state IC) becomes 



l(Oi|C)l' = 



1- 1(^101 



(6) 



which is a linear relation between a square module and a 
module of two probability amplitudes. When the initial 
state is pure, the average of © on the Hilbert space 
reaches the value 3/4. 

Thus, we have found the 9 observable which optimizes 
the fidelity or the probability of taking the initial known 
state p to the target |C) by means of von Neumann mea- 
surements only. 



III. PROCESSING BY JV+1 OBSERVABLES 

First of all we study the case where three observables 
are implemented in order to achieve the target. Hence 
the conlusions obtained are easily generalized when iV + 1 
observables are considered. 
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Now we suppose that before measuring the observable 
C we measure two observables, say 9i followed by 62, 
which define orthonormal bases {|0j), j = 1,2, re- 

spectively. In other words, we shall apply three consec- 
utive von Neumann measurement processes, first M{9i) 
followed by Af(^2) and finally by M{(), which shall be 
denoted by the simple product M(C)M(^2)M(^i). The 
probability of driving the known initial state p toward the 
IC) target, by means of the von Neumann measurement 
M{C)M{e2)M{§i) process, is given by 



P2..S = l-(0i|p|0i)-(l-2(0i|p|0i)) 



i-i(02ic>r-i(0ii02)r i-2|(02ic)i 



(7) 



where the quantities (Oi|p|Oi) and |(0i|02)P entering in 
P2.S are considered to be functions of the coefficients 
of p in the basis of the ^ observable, of the quantities 
|(0i|OP and |(02|C)P, and of the phases ip and of 
(C|p|a)(0i|C)(a|0i) and (0i|O(C|02)(0i|a)(a|02), re- 
spectively. 

The problem of optimizing the probability P2,s, Eq- 
((TJ , leads to a set of nonlinear equations for the quantities 
I (Oi IC) I , I (O2 10 1 if and (j), which can not be analytically 
solved. However, we are able to show that, under certain 
conditions, it is possible to choose the observable 62 in 
such a way that p2,s becomes higher than pi^^- 

The probability Ps^2 can be also written as 



PM + (0i|p|0i)(|(0i|02)n(02|C)|' 
+ |(0i|l2>|'|(l2|C)l'-|(0i|C>l') 
+ (liHli)(|(li|02)n(02|C)|' 

+i(iiii2>ri(i2ic>i'-i(iiic>i'), (8) 

where pi ^ is given by Eq. Hence p2,s is higher than 
Pi,s under the conditions: 



and 



l(0i|02)n(02|C>l' 



l(0i|i2)n(i2|C)l'>l(0ilOl'. 



(9) 



(10) 



The condition I|1U|) means that the basis {IO2), II2)} has 
to be chosen in a way such that the probability of taking 
the state |0i) to the state |C) by means of the M(C)M(6'2) 
process be higher than the probability of taking the state 
|0i) to the state |C) by means of the Af (C) process. We 
have already shown, in section ^ that such a choice is 
always possible. The probability p2.s is higher than pi s 
also under the conditions: 



(OiIpIOi) < (li|p|li) 



(11) 



and 



The latter condition has the same meaning as the 1)10(1 
inequality, but in this case starting from the |li) state 
instead of from the |0i) state. This condition can also be 
always satisfied. 

The above result can be generalized to the case of N 
observables thetai. In this case we suppose that, before 
measuring the observable Ci we measure N observables, 
say 9i, 62, . . ., ^AT, each one defining an orthonormal ba- 
sis respectively, with ij = 0, 1 and j — 1,2 ... ,N. 
The probability of driving the known initial state p to- 
wards the IC) target, by means of the von Neumann mea- 
surement processes M(C)M(^7v) • • ■ M(^2)M(^i), can be 
calculated recursively as 



PN,s = {C\Pn\0: 

where p^ is given by 



(13) 



PN ^ ^ {iN\PN-l\iN)\iN){iN\- (14) 

■ijv=0 

The difference A = p7v+i,s — Pn.s can be read as 



A = ^ {iN\pN-i\iN) 

1 

X ( J2 \{^N\^N+M^N+l\0\'-\{^N\0\')■ 



i]V4-l=0 



(15) 



The positivity of this difference is guaranteed under the 
conditions 



(0Ar|pAr_i|07v) > {1n\PN-i\In) , 



(16) 



and 



|(ll|02)n(02|C)l' + |(ll|l2)n(l2|C)l'>l(ll|Or (12) 



j2 i(o^iz^+i)n(z^+iioi'>i(o^ic)r (17) 

The latter condition means that the basis 
{IOat+i), IIat+i)} must be chosen in a way such 
that the probability of taking state \Qn) to state \() by 
means of the M(C,)M(9n+i) process be higher than the 
probability obtained by means of the M{Q process. In 
section m we have already shown that this can be always 
achieved. The positivity of Eq. (|15|) is also satisfied if 
(OArlpTv-ilOAr) < (1n\pn-i\In) and 

1 

^ \{l^\^^+,)\'\{^^+,\0\' > HUO?, (18) 

which has the same meaning as the H17f) condition, start- 
ing from the IIat+i) state instead of from IOat+i). 

Thus, we have shown that the probability of success 
can be increased by adding suitable observables to the 
process. Since each suitable 9i observable depends on 
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FIG. 3: Maximum value of pjv,s probability as a function of 
the (CIpIC) initial probability for different A'' values: A'^ = 1 
(solid), N = 2 (dash), A'' = 3 (dot), for: (a) an initial pure 
state 7=1 and (b) a mixed initial state 7 = 0. The diagonal 
solid line corresponds to pd- 

the p initial state and the \0 target, we can conclude 
that as N goes to infinity, the fidelity and the probability 
of finding \() will go to 1. This conclusion also can be 
obtained by studying the Hilbert-Schniidt distance 
between pN and |C)(CI- 

Figure |2| shows the results of a numerical simulation 
which finds the probability for a given set of bases 
{|0j>, (j = 1, 2, ... TV) with N and p fixed. On it is 
plotted the maximum value of pn,s as a function of the 
(CIpIC) initial probability for different N values: = 1 
(sofid), A^ = 2 (dash), A^ = 3 (dot), when (a) 7 = 1 and 
(b) 7 = 0. In Fig. |2|a, which corresponds an initial pure 
state, we can see that pn,s increases with respect to pd 
(the diagonal) for all (C|p|C) initial probability. In Fig. 
|31b, which corresponds to an initial mixed state (diagonal 
in the {|C), IC-l)}) basis), we can sec that p^^g increases 
with respect to pd (the diagonal) only for (C|p|C) < 1/2- 

IV. CONCLUSIONS 

In summary, we have studied the problem of driving a 
known initial quantum state onto a known pure state 



without using any unitary transformation. This task 
can be achieved by means of von Neumann measurement 
processes, introducing A'' observables which are consec- 
utively measured in order to get the state closer to the 
target state. We proved that the probability of project- 
ing onto the target can be increased by adding suitable 
observables to the process. Since each of these suitable 
observable depends on the p initial state and on the \0 
target, we conclude that as A^ increase the probability of 
finding \0 goes to 1. 

For a physical implementation of the above described 
process one could address the problem of keeping the ini- 
tial flux of a beam composed of a collection of systems 
in the same state, each one exposed to a postselection- 
measurement procedure. For instance, let us consider a 
source of monochromatic and vertically linear polarized 
photons [23 . In order to obtain photons in a horizontally 
linear polarized state it is required to put a linear polar- 
izer in their path. Implementing two linear polarizers in 
a suitable configuration, the outcome flux with horizontal 
linear polarization is decreased fifty per cent with respect 
to the incoming flux. By implementing more than two 
linear polarizers, as is suggested above, the output flux of 
the beam can be increased meaningfully and it can be ap- 
proached to the initial flux, depending on the number of 
linear polarizers arranged suitably. Since in this scheme 
only one component of the each linear polarized flux con- 
tributes to the success probability, it will converge a little 
more slowly than our protocol, however it will also go to 
1 as the number of linear polarizers arranged suitably 
increases, preserving approximately the initial flux. A 
nonlinear crystal can change the polarization of a pho- 
ton while preserving the flux; however, it also preserves 
the initial mix degree. In our scheme, independently of 
the initial mix degree, the output is pure. 

Further studies could be generalized considering a d- 
dimensional Hilbert space. 

This work was supported by Milenio Grant ICM P02- 
49F, and FONDECYT Grants 1030671, 1040591, and 
1040385. 
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